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Abstract 

The swinging spring, or elastic pendulum, has a 2:1:1 resonance arising at cubic order 
in its approximate Lagrangian. The corresponding modulation equations are the well- 
known three- wave equations that also apply, for example, in laser-matter interaction in 
a cavity. We use Hamiltonian reduction and pattern evocation techniques to derive a 
formula that describes the characteristic feature of this system's dynamics, namely, the 
stepwise precession of its azimuthal angle. 

PACS numbers: 02.40.-k, 05.45.-a, 45.10.Db, 45.20.Jj 

Keywords: Classical mechanics, Variational principles, Averaged Lagrangian, Elastic 
Pendulum, Nonlinear Resonance. 



'email: dholm@lanl.gov 
t email: Peter. Lynch@met.ie 



1 



D. D. Holm & P. Lynch 



Precession of the Swinging Spring 



2 



1 Introduction 

1.1 Problem statement, approach and summary of results 

The elastic pendulum or swinging spring is a simple mechanical system that exhibits 
complex dynamics. It consists of a heavy mass suspended from a fixed point by a light 
spring which can stretch but not bend, moving under gravity. We investigate the 2:1:1 
resonance dynamics of this system in three dimensions and study its characteristic 
feature - the regular stepwise precession of its azimuthal angle. 

When the Lagrangian is approximated to cubic order and averaged over the fast 
dynamics, the resulting modulation equations have three independent constants of 
motion and are completely integrable. These modulation equations are identical to 
the three-wave equations for resonant triad interactions in fluids and plasmas, and 
in laser-matter interaction. We reduce the system to a form amenable to analytical 
solution and show how the full solution may be reconstructed. We examine the 
geometry of the solutions in phase-space and develop a number of simple qualitative 
descriptions of the motion. 

We compare solutions of the exact and modulation equations and show that they 
are remarkably similar. A characteristic stepwise precession occurs as the motion 
cycles between quasi-vertical and quasi-horizontal. That is, during each quasi-vertical 
phase, the azimuth of the swing plane precesses by a constant angular increment. This 
stepwise azimuthal precession occurs in bursts, when the motion is nearly vertical. 
By transforming to non-uniformly rotating coordinates and assuming a geometric 
constraint (essentially the method of pattern evocation), we find a formula for the 
rotation of the swing plane. This formula gives a highly accurate description of the 
stepwise precession of the azimuthal angle of the motion. 

1.2 History of the problem 

The first comprehensive analysis of the elastic pendulum appeared in Vitt and Gorelik 
(1933). These authors were inspired by the analogy between this system and the Fermi 
resonance of a carbon-dioxide molecule. We make connections in this paper with other 
physical systems of current interest. For example, we show that the modulation 
equations for the averaged motion of the swinging spring may be transformed into 
the equations for three-wave interactions that appear in analyzing fluid and plasma 
systems, and in laser-matter interaction. These three complex equations are also 
identical to the Maxwell-Schrodinger envelope equations for the interaction between 
radiation and a two-level resonant medium in a microwave cavity (Holm and Kovacic, 
1992). The three- wave equations also govern the envelope dynamics of light waves in 
an inhomogeneous material (David et al, 1990, Alber et al, 1998a,b). For the special 
case where the Hamiltonian takes the value zero, the equations reduce to Euler's 
equations for a freely rotating rigid body. Finally, the equations are also equivalent to 
a complex (unforced and undamped) version of the Lorenz (1963) three-component 
model, which has been the subject of many studies (Sparrow, 1982). Thus, the 
simple spring pendulum, which was first studied to provide a classical analogue to 
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the quantum phenomenon of Fermi resonance, now provides a concrete mechanical 
system which simulates a wide range of physical phenomena. 

All of the previous studies of the spring pendulum known to us have considered 
motion in two dimensions. To our knowledge, only Cayton (1977) discussed three- 
dimensional solutions and observed the curious rotation of the swing plane between 
successive cycles when the horizontal energy is maximum. This particular aspect of 
the behavior of the swinging spring in three dimensions is its most striking difference 
from two-dimensional motions. Suppose the system is excited initially near its vertical 
oscillation mode. Since purely vertical motion is unstable, horizontal motion soon 
develops. The horizontal oscillations grow to a maximum and then subside again. An 
alternating cycle of quasi-vertical and quasi-horizontal oscillations recurs indefinitely. 
Seen from above, during each horizontal excursion of several oscillations the projected 
motion is approximately elliptical. Experimentally and numerically one observes that 
between any two successive horizontal excursions the orientation of the projected 
ellipse rotates by the same angle, thereby causing a stepwise precession of the swing 
plane. In principle, the precession angle between successive horizontal excursions can 
be deduced from the complete solution of the integrable envelope equations. We seek 
a simple approximate expression for the precession of the swing plane in terms of the 
solution of the reduced dynamics. 

Lynch (2001) found a particular solution for the rate of precession of the swing 
plane by using the method of multiple time scales in rotating coordinates and in- 
troducing a certain angular solution Ansatz. We recover Lynch's particular solution 
among a family of other solutions for the swing plane precession rate. This fam- 
ily is obtained via the method of averaged Lagrangians by seeking solutions of the 
modulation equations that satisfy a geometrical constraint of being "instantaneously 
elliptical." We apply the method of pattern evocation in shape space (Marsden, et 
al, 1995, 1996). Using this process, one identifies patterns by viewing the dynamics 
relative to rotating frames with certain critical angular velocities. Our numerical in- 
tegrations show that the solution resulting from this geometrical postulate estimates 
the precession of the swing plane with surprisingly high accuracy. 

2 Equations of motion 

The physical system under investigation is an elastic pendulum, or swinging spring, 
consisting of a heavy mass suspended from a fixed point by a light spring which can 
stretch but not bend, moving under gravity, g. We assume an unstretched length £ , 
length £ at equilibrium, spring constant k and unit mass m — 1. The corresponding 
Lagrangian, approximated to cubic order in the amplitudes, is 

L = \ [x 2 + y 2 + z 2 ) - \ (u 2 R (x 2 + y 2 ) + cu 2 z 2 ) + \\{x 2 + y 2 )z , (2.1) 

where x, y and z are Cartesian coordinates centered at the point of equilibrium, 
ojr = \fg~jt~ is the frequency of linear pendular motion, u>z = \Jkjm is the frequency 
of its elastic oscillations and A = ^uj 2 z /£ 2 . The system is illustrated schematically in 
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Fig. 1. The Euler-Lagrange equations of motion may be written 



x + UJ R X 

V + u 2 R y 

Z + ijJ 2 7 Z 



Xxz , 
Xyz, 



(2.2) 
(2.3) 
(2.4) 



There are two constants of the motion, the total energy E and the angular momentum 
h given by 

E = \{x 2 + y 2 + z 2 ) + \{u R {x 2 + y 2 ) + to 2 z z 2 ) - \X{x 2 + y 2 )z , h = {xy - yx) . 

The system is not integrable. Its chaotic motions have been studied by many authors 
(see, e.g., Refs. in Lynch, 2004). Previous studies have considered the two dimensional 
case, for which the angular momentum h vanishes. 

We confine attention to the resonant case uz = and apply the averaged 
Lagrangian technique (Whitham, 4974). The solution of fl2.2|) -( f2~4"D is assumed to be 
of the form 



x = 3?[ao(i) exptyujRt)] , 
y = R[b (t)exp(iuj R t)], 
z = 3?[c (t) exp(2iuj R t)} . 



(2.5) 
(2.6) 
(2.7) 



(The zero-subscripts in a , b and c are introduced to distinguish from the variables 
a, b and c in a rotating frame, introduced below.) The coefficients a (t), bo(t) and 
co(t) are assumed to vary on a time scale which is much longer than the time-scale of 
the oscillations, r = 1/uj r . The Lagrangian is averaged over time r to give, 



(L) 



5^R 



Q{a a* + b b* + 2c c* } + K$t{(a 2 + b 2 )c* } 



where k = X/(4uj r ). We regard the quantities a ,&0) c o as generalized coordinates. 
The averaged Lagrangian equations of motion are then 



ib 
ic 



na Q CQ , 
Kb* c , 



(2.8) 
(2.9) 
(2.40) 



Eqns. ( |2.8p -( |2~TUD are the complex versions of Eqns. (68)-(73) in Lynch (2004), which 
were derived using the method of multiple time-scale analysis. 

We now absorb the constant k by rescaling the time, t — > Kt, and transform 
variables 



A = ±(a + ib ), B = ±(a -ib ), 

Consequently, the equations of motion take the form 

iA = B*C, 
iB = CA\ 
iC = AB. 



c 



c - 



(2.11) 
(2.12) 
(2.13) 
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These three complex equations are well-known as the three-wave interaction equations, 
which govern quadratic wave resonance in fluids and plasmas. 

The three-wave interaction equations ( 2.11|) -( p.l3|) may be written in canonical 



form with Hamiltonian H = $l(ABC*) and Poisson brackets {A, A*} = {B,B*} 
{C,C*} = -2i, as 

iA = i{A, H } = 2dH/dA* , (2.14) 
iB = i{B, H} = 2dH/dB* , (2.15) 
iC = i{C, H} = 2dH/8C* . (2.16) 

These equations conserve the following three quantities, 

H = \{ABC* + A*B*C) = $1{ABC*) , (2.17) 
|2 + 2|Cf, (2.18) 

(2.19) 



N = 


\A\ 2 + 


\B\ 


J = 


\A\ 2 - 


\B\ 
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Thus, the modulation equations for the swinging spring are transformed into the 
three-wave equations, which are known to be completely integrable. See Alber et al. 
(1998a) for references to the three wave equations and an extensive elaboration of 
their properties as a paradigm for Hamiltonian reduction. 

The following positive-definite combinations of iV and J are physically significant: 

N + = i(iV + J) = \A\ 2 + \C\ 2 , iV_ = l(N -J) = \B\ 2 + |Cf . 

These combinations are known as the Manley-Rowe relations in the extensive liter- 
ature about three-wave interactions. The quantities H, N + and AL provide three 
independent constants of the motion. 

2.1 A brief history of the three- wave equations 

Fluids and plasmas. The three-wave equations model the nonlinear dynamics of 
the amplitudes of three waves in fluids or plasmas (Bretherton, 1964). The equations 
result from a perturbation analysis of the Charney equation 

& tTj2 i 7n;i , fdtpdV 2 ^ dtp dV 2 ip\ t difj 



^ [V ^ ^ ydx dy dy dx ) ~*~ ^ dx ^ ^ ^ 

(see, e.g., Pedlosky, 1987 for theory and notation). This equation is equivalent to 
the Hasegawa-Mima equation describing drift-waves in an inhomogeneous plasma 
in a magnetic field (Hasegawa and Mima, 1977). Longuet-Higgins and Gill (1967) 
examined the interactions between planetary Rossby waves in the atmosphere and 
derived detailed conditions for three-wave resonance. The correspondence between 
Rossby waves in the atmosphere and drift waves in plasma have been thoroughly 
explored by Horton and Hasegawa (1994). Resonant wave-triad interactions play an 
essential role in the generation of turbulence and in determining the statistics of the 
power spectrum. Both energy and enstrophy are conserved in fluid systems governed 
by the Charney equation ( |2.20| ). 
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Laser-matter interaction. Equations ( j2.11 )-( 2.13| ) are also equivalent to the 



Maxwell-Schrodinger envelope equations for the interaction between radiation and 
a two-level resonant medium in a microwave cavity. Holm and Kovacic (1992) show 
that perturbations of this system lead to homoclinic chaos, but we shall not explore 
that issue here. Wersinger et al. (1980) used a forced and damped version of the 
three-wave equations to study instability saturation by nonlinear mode coupling, and 
found irregular solutions indicating the presence of a strange attractor. See also Ott 
(1993) and Holm, Kovacic and Wettergren (1995, 1996) for more detailed studies of 
the perturbed three-wave system. 



Nonlinear optics. The three-wave system also describes the dynamics of the en- 
velopes of light-waves interacting quadratically in nonlinear material. The system 
has been examined in a series of recent papers (Alber, et al, 1998a,b; Luther et al., 
2000) using a geometrical approach which allowed the reduced dynamics for the wave 
intensities to be represented as motion on a closed surface in three dimensions - the 
three-wave surface. Information about the corresponding reconstruction phases was 
recovered using the theory of connections on principal bundles. 

In the special case H = the system fl2.11| )-( p.l3p reduces to three real equations. 

Let 

A = iXiexp^x) , B = iX 2 exp(i(j) 2 ) , C = zX 3 exp(i(0x + 2 )) 

where X\, X 2 and X 3 are real and the phases and cj) 2 are constants. The modulation 
equations become 

X\ = — X 2 X 3 , X 2 = — X 3 X 1 , X 3 = +XiX 2 . (2-21) 

We note that these equations are re-scaled versions of the Euler equations for the 
rotation of a free rigid body. The dynamics in this special case is expressible as 
motion on R 3 namely, 

X = f V J x ViV = ±ViV + x ViV_ . (2.22) 

Considering the constancy of J and N, we can describe a trajectory of the motion as 
an intersection between a hyperbolic cylinder (J constant, see (|2.19| )) and an oblate 
spheroid (N constant; see (|2.18|) ). Eqn. ( |2.22|) provides an alternative description. 
Here we have used the freedom in the R 3 Poisson bracket exploited by David and 
Holm (1992) to represent the equations of motion on the intersection of two orthogonal 
circular cylinders, the level surfaces of the Manley-Rowe quantities, iV + and iV_. 
The invariance of the trajectories means that while the level surfaces of J and iV 
differ from those of iV + and AL, their intersections are precisely the same. For this 
particular value of H = 0, the motion may be further reduced by expressing it in the 
coordinates lying on one of these two cylinders. See Holm and Marsden (1991) for the 
corresponding transformation of rigid body motion into pendular motion. See David 
and Holm (1992), and Alber et al. (1998a,b) for discussions of geometric phases in 
this situation. 
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2.2 Reduction of the system & Reconstruction of the solution 

To reduce the system for H 7^ 0, we employ a further canonical transformation, 
introduced in Holm and Kovacic, 1992. Namely, we set 

A = |A|exp(if), (2.23) 
B = |5|expO?), (2.24) 
C = Z exrj(i(£ + ri)) . (2.25) 

This transformation is canonical - it preserves the symplectic form 

dA A dA* + dB A dB* + dC A dC* = dZ A dZ* . 

In these variables, the Hamiltonian is a function of only Z and Z* 

H = \{Z + Z*) ■ v/iV+ - \Z\ 2 ■ y/N_ - \Z\ 2 . 

The Poisson bracket is {Z, Z*} = —2i and the canonical equations reduce to 

rll-f 

iZ = i{Z,H} = 2—. 

This provides the slow dynamics of both the amplitude and phase of Z = \Z\e 1 ^. 

The amplitude \Z\ = \C\ is obtained in closed form in terms of Jacobi elliptic 
functions as the solution of 

(It) = [2 3 - 2 2 2 + ( 1 -^ 2 )2 + 2 ^]' ( 2 - 26 ) 

where Q = 2\Z\ 2 /N, J = J/N, £ = —AH 2 /N 3 and r = y/2Nt. This is equivalent to 
Eqn. (75) in Lynch, 2001 (an explicit expression for the solution in terms of elliptic 
functions is given in that paper). Once \Z\ is known, \A\ and \B\ follow immediately 
from the Manley-Rowe relations, 



\A\ = - \Z\ 2 , \B\ = v/iV_ - \Z\ 2 . 

The phases ^ and r\ may now be determined. Using the three-wave equations ( |2.11 )- 
( glg ) together with (ggg ) - (ggg) , one finds 

i = ~W> i, = ~W' (2,27) 

so that ^ and rj can be obtained by quadratures. Finally, the phase ( of Z is deter- 
mined unambiguously by 

2 

= -2#tanC and H = \A\\B\\Z\ cos( . (2.28) 



d\Z\ 



dt 

Hence, we can now reconstruct the full solution as 



A = \A\ exp(i^) , B = \B\exp(ir)) , C = \Z\ exp (z(£ + r/ + ()) . 



D. D. Holm & P. Lynch 



Precession of the Swinging Spring 



8 



3 Phase portraits. 

Consider the plane C in phase-space defined by A = B = 0. This is a plane of 
unstable equilibrium points, representing purely vertical oscillations of the spring. 
The Hamiltonian vanishes identically on this plane, as does the angular momentum 
J. Each point Cq in C has a heteroclinic orbit linking it to its antipodal point — cq. 
Thus, the plane C of critical points is connected to itself by heteroclinic orbits. In 
Fig. 2, the horizontal plane is C and the vertical plane contains heteroclinic orbits 
from Co to — cq. The vertical axis is R = a/I^I 2 + \B\ 2 . Since N = R 2 + 2\C\ 2 is 
constant, each heteroclinic orbit is a semi-ellipse. Motion starting on one of these 
semi-ellipses will move towards an end-point, taking infinite time to reach it. 

In Fig. 3 taken from Holm and Kovacic 1992, we present another view of the 
trajectories for J = 0. The Hamiltonian is 

H=\{Z + Z*) ■ {\N- \Z\ 2 ). 

Accessible points lie on or within the circle \Z\ 2 = N/2. For H = the trajectory is 
the segment of the imaginary axis within the circle. This is the homoclinic orbit. For 
H 0, we solve for the imaginary part of Z = Z\ + zZ 2 , 

Z 2 = ±sJ-Z 2 + \N -{H/ZJ 

This allows us to plot the trajectories for the range of H for which real solutions exist. 
There are two equilibrium points, at Z = ±a/ N/6, corresponding to solutions for 
which there is no exchange of energy between the vertical and horizontal components. 
These are the cup-like and cap-like solutions first discussed by Vitt and Gorelik (1933). 



3.1 Geometry of the motion for fixed J 

The vertical amplitude is governed by equation ( |2.26| ), which we write as 



1 fdQ\ 



2 



2{W) +V(S)=£ ' M 
with the potential V(Q) given by 

V(Q) = -\ [Q 3 - 2Q 2 + (1 - J 2 )Q] . (3.2) 



We note that V(Q) has three zeros, Q = 0, Q = 1 - J and Q = 1 + J. Eqn. Q is 
an energy equation for a particle of unit mass, with position Q and energy £, moving 
in a cubic potential field V(Q). In Fig. 4 we plot Q, given by (ft.l|) , against Q for 
the cases J = (left panel) and J = 0.25 (right panel), for a range of values of £. 
Each curve represents the projection onto the reduced phase-space of the trajectory 
of the modulation envelope. The centers are relative equilibria, corresponding to the 
elliptic-parabolic solutions of (Lynch, 2001), which are generalizations of the cup- 
like and cap- like solutions of Vitt and Gorelik (1933). The case J = includes the 
homoclinic trajectory, for which H = 0. 
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3.2 Geometry of the motion for H = 

For arbitrary J , the H = motions are on a surface in the space with coordinates 
(Q, Q, J)- This surface is depicted in Fig. 5. It has three singular points (i.e., it is 
equivalent to a sphere with three pinches) and its shape is similar to a tricorn hat. 
The motion takes place on an intersections of this surface with a plane of constant J . 
There are three equilibrium solutions: that with J = (marked H.P. in Fig. 5) is at 
the extremity of the homoclinic trajectory, and corresponds to purely vertical oscil- 
latory motion; those with J = ±1 correspond to purely horizontal motion, clockwise 
or anti-clockwise, with the spring tracing out a cone. The purely vertical motion is 
unstable; the conical motions are stable (perturbations about conical motion were 
investigated by Lynch, 2001). The dynamics on the tricorn are similar to the motion 
of a free rigid body. The three singular points correspond to the steady states of 
rotation about the three principal axes. 



3.3 Three- wave surfaces 

There is another way to depict the motion in reduced phase-space. Let us consider 
a reduced phase-space with x and y axes X = dt{ABC*} and Y = Q{ABC*} and 
z-axis Q = 2\Z\ 2 /N. We note that X = H. It follows from (PT7Q-(P^|) that 



X 2 + Y 2 = \A\ 2 \B\ 2 \C\ 2 = \\Z 



{2\Z\ 2 — N) 2 — J 2 

We define X = (2/N 3 / 2 )X and y = (2/N 3 l 2 )Y and can write 

X 2 + y 2 = \ [Q 3 - 2Q 2 + (1 - J 2 )Q] = -V(Q) . (3.3) 

where V is defined in We note that X 2 = -8 and y 2 = ^(dQ/dr) 2 . Eqn. Q 

implies that the motion takes place on a surface of revolution about the Q-axis. The 
radius for a given value of Q is the square-root of the cubic — V(Q). The physically 
assessable region is < Q < 1 — \J~\. Several such surfaces (for J G {0.0, 0.1, 0.2, 0.3}) 
are shown in Fig. 6. Since X 2 = 7i 2 = 4H 2 /N 3 , the motion for given J takes place 
on the intersection of the corresponding surface of revolution with a plane of constant 
X. 

We can relate the tricorn surface to the surface of revolution. The former is 
appropriate for if = 0; the H ^ case is represented by trajectories inside this 
surface. If we slice the tricorn surface in a plane of fixed J we get a set of closed 
trajectories, the outside one for H = the others for H ^ (the cases J = and 
J = 0.25 are plotted in Fig. 4 above.) If we now distort the j7"-section into a cup-like 
surface, by taking H as a vertical coordinate and plotting each trajectory at a height 
depending on its H value, we get half of a closed surface. Each trajectory is selected 
by a if-plane section. Alteration of the sign of H corresponds to reversal of time. 
Completing the surface by reflection in the plane H = gives the surface generated 
by rotating the root-cubic graph y — V(Q) about the Q-axis, i.e., the surface given 
by (|3.3|) . These surfaces are what Alber, et al, (1998a) call the three- wave surfaces. 
They foliate the volume contained within the surface for J = 0. 
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4 The Precession of the swing plane 

The characteristic feature of the behavior of the physical spring is its stepwise preces- 
sion, which we shall now analyze. As the oscillations change from horizontal to vertical 
and back again, it is observed that each successive horizontal excursion departs in a 
different direction. The only reference to this phenomenon of which we are aware, 
prior to Lynch (2001), is Cayton (1977). Cayton briefly discussed three-dimensional 
solutions and mentioned the precession of the swing plane, but did not analyze its dy- 
namics. Surprisingly, the characteristic stepwise precession of the swinging spring has 
been largely ignored, although it is immediately obvious upon observation of a phys- 
ical elastic pendulum with uj z 2uj r . Indeed, this precession is almost impossible to 
suppress experimentally when the initial motion is close to vertical. 

4.1 Qualitative description 

If the horizontal projection of the motion is an ellipse of high eccentricity, the motion is 
approximately planar. We call the vertical plane through the major axis of this ellipse 
the swing plane. When the initial oscillations are quasi-vertical, the motion gradually 
develops into an essentially horizontal swinging motion. This horizontal swinging does 
not persist, but soon passes again into nearly vertical springing oscillations similar 
to the initial motion. Subsequently, a horizontal swing again develops, but now in 
a different direction. The stepwise precession of this exchange between springing 
and swinging motion continues indefinitely in the absence of dissipation and is the 
characteristic experimental feature of the swinging spring. We shall seek an expression 
for the change in direction of the swing plane from one horizontal excursion to the 
next. 

A full knowledge of the solutions of the three equations of motion would of course 
suffice to determine the swing plane at each moment in time. In Lynch (2001) the 
equations were expressed in rotating co-ordinates, and a particular solution for the 
slow rotation of the swing plane was posited as a function of the vertical amplitude 
\C\, by assuming a certain angular relation. Following this assumption, the angle of 
the swing-plane could be expressed as an integral involving elliptic functions. 

4.2 Pattern evocation in shape space 

We shall approach the precession problem using pattern evocation in shape space. 
Pattern evocation seeks a relative equilibrium (in shape space) in which a phase 
relationship between the variables (the shape) is preserved (Marsden, et al, 1995, 
1996). We track the pattern by moving to a non- uniformly rotating frame in which 
the orientation of the shape is fixed. This is a generalization of the idea of tracking a 
satellite orbit by evoking constancy of the areal velocity required to conserve angular 
momentum. 

Our particular geometric assumption is that the angle between the complex am- 
plitudes a and b remains constant, in an appropriately rotating frame. Writing these 
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amplitudes in vector form as a = (\a\ cos a, \a\ sina,0), b = (|6| cos/3, |6| sin/5,0) and 
taking k = (0, 0, 1) yields 

J = — k ■ a x b = \ab\ sin(o; — j3) , a ■ b = \ab\ cos(a — (3) . 

Consequently, our geometric pattern evocation assumption that the phase difference 
a — (3 remains constant immediately implies that \ab\ is also constant. The conser- 
vation of angular momentum J means that the area of the parallelogram formed by 
the vectors a and b is constant. The requirement of constant a — (3 imposes an ad- 
ditional geometric constraint on the possible shape of the orbits. For example, when 
a — (3 = 7r/2 (mod7r), the orbits are elliptical. 

4.3 Modulation equations in rotating coordinates 

We shall transform to rotating co-ordinates and seek an expression for the (slow) 
rotation frequency Q(t) that allows us to estimate the stepwise precession of the 
swinging spring by imposing the pattern evocation constraint that a — (3 remains 
constant. 

In a rotating frame, the approximate Lagrangian ( |2.1|) at cubic order in the 
coordinate displacements becomes, with x = (x,y,z), 

L = i|x + Q(t) z x x| 2 - | (c4(x 2 + y 2 ) + cu 2 z z 2 ) + \\{x 2 + y 2 )z . (4.1) 

Now x, y and z are Cartesian coordinates centered at the point of equilibrium in the 
rotating frame, lor = \J gjl is the frequency of linear pendular motion, uz = \fkfm 
is the frequency of its elastic oscillations and A = iou^/i 2 . The corresponding Euler- 
Lagrange equations of motion ( |2.2|) -( pT4|) may be written in rotating coordinates as 

x - tl(t)y - 2Q(t)y + (oo 2 R - Q 2 (t)) x = Xxz , (4.2) 
y + Q{t)x + 2Q{t)x + (u 2 R -Q 2 (t))y = Xyz , (4.3) 

z + u 2 z z = \\{x 2 + y 2 ). (4.4) 

The vertical component of angular momentum is 

h = z • x x (± + Q(t) z x xj = (xy - xy) + Q(t)(x 2 + y 2 ) 

and is a constant of the motion for these equations. However, upon Legendre- 
transforming, one finds the time-dependent Hamiltonian satisfies 

H = -tt(t)h. 

Thus, perhaps not unexpectedly, exact energy conservation breaks down to the extent 
that the rotation frequency is nonuniform. 
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4.4 Averaged Lagrangian and modulation equations for slow rotation 



The modulation equations in axes rotating with angular velocity Q(t) about the verti- 
cal are obtained in the resonant case u>z = 2lv r by applying the averaged Lagrangian 
technique (Whitham, 1974). Accordingly, the solution of ( f4.2|) -( f4~4"|) is assumed to be 
of the form 



x 

y 

z 



fft[b(t)esxp(iu R t)], 
U[c(t) exp(2icj R t)} . 



(4.5) 
(4.6) 
(4.7) 



(Note that subscript-zeroes are dropped for these modulation amplitudes in the ro- 
tating frame.) In these variables, the averaged Lagrangian corresponding to 
may be written as 



(L) 



<3{aa* + bb* + 2cc*} + 3?{(a 2 + b 2 )c*} + 2Q%{ab*} 
+ if) 3? [a*b - a*b] + \tt 2 \\a\ 2 + \b\ 



(4i 



On assuming that the rotation frequency is sufficiently slow that VL/ujr <C 1, we shall 
neglect all terms in the averaged Lagrangian ([4.8| ) that are not multiplied by uor. In 
this approximation of slow rotation, the averaged Lagrangian is given by the simpler 
expression, 



(L) 



%{aa* + bb* + 2cc*} + 9tf{(a 2 + b 2 )c*} + 2QJ 



(4.9) 



Here J = Q{ab*} is the angular momentum, a conserved quantity at this level of 
approximation and formally identical to the expression in non-rotating coordinates. 
The Euler-Lagrange modulation equations in this approximation may be written as 



ia 
ib 
it 



a*c + iVtb . 
b*c — ifla . 



+ b 2 ). 



(4.10) 
(4.11) 
(4.12) 



We may also write these leading order equations in Hamiltonian form. When (H) is 
defined by 

(H) = 3?{(a 2 + b 2 )c*} + 2QQ{ab*} , 

with coordinates (a, b, c), conjugate momenta (a*, b*, 2c*) and Poisson brackets defined 
by {a, a*} = {b,b*} = 2{c, c*} = —i, the modulation equations ( |4.10p -( P~T2"D are 
expressible in canonical Hamiltonian form as, 



i a 



i{a,(H)} 



d(H) 
da* 



ib = i{b,(H)} 



9(H) 
db* 



ic 



z{c, (H)} 



d(H) 
02c* 
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The three constants of the motion for these equations are 

H = l[( a 2 + b 2 )c* + (a* 2 + b* 2 )c} = ^{(a 2 + b 2 )c*}, (4.13) 
N = \a\ 2 + \b\ 2 + 4|c| 2 , (4.14) 
J = (ab* -a*b)/2i = S{ab*} . (4.15) 

We now introduce the pattern evocation assumption. Noting that 

\ab\ 2 = (^{ab*}) 2 + (%{ab*})\ 

and that the second term is just J 2 , implies constancy of 3?{a6*}. Using Eqns. Q4.10|) 
and (|4.11| ) it follows that 

j t \ab\ 2 = -2^{ab*} [2%{abc*} + ft(|a| 2 - |6| 2 )] = . (4.16) 

Either factor may vanish, so there appear to be two possibilities for the solution. We 
first assume that the factor in square brackets in (|4.16|) vanishes. This implies 

_ 2%{abc*} _ 2\abc\ sm(a + g - 7 ) 

\a\ 2 -\b\ 2 ~ |a| 2 -|6| 2 ' [ ' 

(where c = |c|e 17 ). The precession ang le 6 = j Q VL{t')dt' can be ascertained by 
integrating Q over the time interval of the motion. In the special case a — (3 = 
|(mod7r), one finds by using the constants of motion that 

"= (Jy -4|ffi-4/. - (4 - 18) 

This case also corresponds to the vanishing of the first factor in (|4.16|) , so that 
3l{ab*} = and a and b are 90° out of phase. This was the Ansatz introduced 
by Lynch (2001). He showed that, in this case, the rotation rate is given by (|4.18| ). 
We now see that the result in Lynch (2001) is a special case of the general result 
( (4.17p . m tins special case, Q can be computed as soon as |c| is known. We will 
examine this case numerically below. 



4.5 The instantaneous ellipse 

In order to define precisely the precession angle, we introduce an ellipse which ap- 
proximates the horizontal projection of the trajectory of the pendulum. Recall that 
the full solution for the horizontal components is 

x = ?R{aexrj(iu!fit)} = \a\ cos(u R t + a) , y = ?R,{bexp(iujjit)} = \b\ cos(o;#t + (3) , 

where a and j3 are the phases of a and b. The amplitudes and phases are assumed 
to vary slowly. If they are regarded as constant over a period r = 1/ujr of the fast 
motion, these equations describe a central ellipse, 



Px 2 + 2Qxy + Ry 2 = S 



(4.19) 
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where P = \b\ 2 , Q = —\ab\cos(a — 8), R = \a\ 2 and S = J 2 . The area of the 
ellipse is easily calculated and is found to have the constant value ir J. Its orientation 
is determined by eliminating the cross-term in ( f4.19| ). This is achieved as usual by 
rotating the axes through an angle 9, given by 

2Q 2\ab\ cosfct — 8) 

***> = p=h= H.-| t p • < 4 - 20 > 

The semi-axes of the ellipse are given by 

J * J 

Ai = ===== ===== , A 



y/Pcos 2 9 + Qsin29 + Rsm 2 9 ' " y/P sin 2 9 - Q sin 29 + R cos 2 9 ' 

(4.21) 

The area is ttAiA 2 = ttJ and the eccentricity can be calculated immediately. In the 
case of unmodulated motion, such as the elliptic-parabolic modes described in Lynch 
(2001), the instantaneous ellipse corresponds to the trajectory, which is a precessing 
ellipse. In general, it is only an approximation to the trajectory, but we may define 
the orientation or azimuth at any time to be the angle 9 given by ( [4.20 ). This angle 



will be compared to the precession angle calculated by integrating Q4.18 ) and shown 
to give almost identical results. 



5 Numerical results 

We examine the results of numerical integrations of the modulation equations 
( |2.1U| ) and compare them to the solutions of the exact equations (|2.2])-(|2~^|). It will be 



seen that the modulation equations provide an excellent description of the envelope 
of the rapidly varying solution of the full equations. We then compare the stepwise 
precession angle predicted by a formula based on constancy of the angle a — 8 with 
the numerical simulation of this quantity and show that the two values track each 
other essentially exactly. 

The parameter values chosen for all numerical integrations are m = 1 kg, i = 1 m, 
g = 7r 2 m s~ 2 and k = 4it 2 kg s~ 2 so that ojr = 7T, ujz = 2tt and the resonance condition 
ujz = 2ujr holds. The linear rotational mode has period = 2 s, and the vertical 
mode has period t% — 1 s. The initial conditions are set as follows, 

(x , y , z ) = (0.006, 0, 0.012) ; (x , y , z ) = (0, 0.00489, 0) . 

(The value of yo was chosen to tune the precession angle to be an even fraction of 
180°, making the amplitudes, though not the phases, periodic). The corresponding 
initial values for the modulation equations (|2.8|) - ([2.1Q ) are given by 



arctan^V & = arctan V 7o = arctan ( 

\urx J \URyoJ \2ujrZqJ 



V 2u)rZq J 

cosa o y ' ' Ul \ursoi 8q ) 1 1 Ul \yCOS7o, 



«o| = , N = - — r r > co 
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giving the values (|a |, |6 |, |c |) = (0.006,0.002,0.012) and (a ,A),7o) = (0, -vr/2,0). 
The constants of the motion take the following values, 

H = 4.03 x 10" 7 , J = 9.34xl0" 6 , N = 6.14 x 10" 4 . 

The integration was extended over a period of 1000 seconds (i.e., 1000 vertical oscil- 
lations). As a check on numerical accuracy, the changes in these quantities, which 
should remain constant, was calculated, with the following results: 

(^ L ) = 100.04% , (^™L) = 99.997% , ( = 100.00% . 

V " Initial / V -'Initial / \ ^Initial / 



We now directly compare the solutions of the 'exact' equations fl2.2|) - (|2.4j) and the 
'approximate' or modulation equations ( ^~8|) -( ^~l0| ). Once the modulation equations 
have been solved for the envelope amplitudes and phases, the full approximate solu- 
tion is given by (|2.5|) - (|2.7|) . We first consider the horizontal projection of the solution 
for the 1000 second integration. This is the period required for the solution to precess 
through approximately 180°. In Fig. 7 (top panel) we plot x versus y for the exact 
solution. In Fig. 7 (bottom panel) we plot the corresponding solution from the mod- 
ulation equations. It is clear that there is great similarity between the two solutions; 
indeed, the two plots are indistinguishable. The precession angle between horizontal 
excursions is close to 30° (the value of yo was chosen to ensure this). The modulation 
period is approximately 167 seconds; thus, the instantaneous ellipse rotates through 
six cycles and 180° in 1000 seconds. 

The vertical structure of the solution is displayed in Fig. 8, where z for the exact 
solution (top panel) and $t{c (t) exp(2iuj R t)} for the approximate solution (bottom 
panel) are seen to be virtually identical. For clarity, the solutions are plotted only 
for the first modulation cycle of 167 seconds. The character of the solution — rapid 
oscillations with a slowly-varying amplitude envelope - - is clear from the figure. 
The vertical amplitude is close to zero when horizontal excursions are at their peak. 
This is confirmed in Fig. 9 (top panel) where the horizontal modulation amplitude 
S = \/\a\ 2 + \b\ 2 and vertical modulation amplitude C — |c| are plotted against time. 



In Fig. 9 (bottom panel) we plot the squared eccentricity e 2 = (1 — A\ A J A 2 



maj/ 



of the envelope of the horizontal projection of the approximate solution, where the 
semi-axes A ma j and A min are calculated from fl4.21|) . The eccentricity is close to unity 



for most of the integration. Horizontal excursions of the pendulum occur during this 
time. For short periods, when the horizontal amplitude is minimum, the value of e 
drops significantly (solid line). During this time, the angular velocity, calculated as 
the rate of change of the azimuth given by (|4.20|) , reaches a maximum (dashed line). 
Thus, the precession occurs in bursts near the times when the vertical amplitude is 
maximum and horizontal amplitude minimum. 

The stepwise nature of the precession is clearly illustrated in Fig. 10. The az- 
imuthal angle $ of the numerical solution of the exact equations may be calculated 
by fitting a central conic to every three consecutive points on the trajectory. Assuming 
a solution of the form 

Px 2 + 2Qxy + Ry 2 = 1 
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and requiring that the_ three points lie on this curve, we obtain three equations for 
the coefficients (P, Q, R). ^From these, the azimuth $ and the semi-axes are obtained 
from equations analogous to ( |4.20| ) and ( |4.21| ). This is compared in Fig. 10 to the 
corresponding value 9 resulting from integration of the modulation equations. It is 
noteworthy that i? and 9 remain quasi-constant for most of the modulation cycle, 
changing rapidly only over short intervals around the times when C is maximum and 
S is minimum. The advances in phase are very similar for the exact and approximate 
solutions. However, there are small differences: 9 — $ is also plotted in Fig. 10 (dotted 
line). This sensitive quantity reaches its maximum value of 4.35° at the end of the 
integration. 

The azimuthal change between successive horizontal excursions is very close to 30° 
for both exact and approximate solutions. We also calculated the angle O resulting 
from an integration of (|4.18| ). The graphs of 9 and (not plotted) are indistinguish- 
able. The maximum difference \9 — 0| was only 0.0063°. This is remarkable: the 
value derived from ([4.18 ) involves an assumption that a — (3 is constant in a par- 
ticular rotating frame. The azimuth 9 from the modulation equations makes no such 
assumption, yet the two solutions are practically identical. This confirms that the 
pattern evocation assumption which yields the result ( |4. 18| ) is sound. 

Numerous other integrations of the exact and modulation equations were also 
carried out. They confirm that the stepwise precession of the azimuthal angle is a 
distinct characteristic of the swinging spring. This is also in complete agreement 
with simple experiments with a physical pendulum, where the periodic exchange of 
energy between horizontal and vertical and the precession of the swing plane between 
horizontal excursions are the main observable properties of the motion. 
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Figure Captions 

Figure 1. Schematic diagram of the elastic pendulum or swinging spring. Cartesian 
coordinates centered at the position of equilibrium are used. 

Figure 2. C is the plane of critical points, A = = B. The vertical axis is R = 
^/\A\ 2 + \B\ 2 . The vertical plane contains heteroclinic semi-ellipses passing from 
c to -Co. 

Figure 3. Phase portrait in the Z-plane for J = 0. The motion is confined within 
the circle \Z\ 2 = ^N. The segment of the imaginary axis within this circle is the 
homoclinic orbit. 

Figure 4. Plots of Q versus Q for J = and J = 0.25, for a range of values 
E e {-0.0635, -0.0529, -0.0423, -0.0317, -0.0212, -0.0106, 0}. 

Figure 5. Tricorn surface, upon which motion takes place when H = 0. The coor- 
dinates are J y Q, Q. The motion takes place on the intersections of this surface 
with a plane of constant J (such planes are indicated by the stripes). This surface 
has three singular points. The homoclinic point is marked H.P. 

Figure 6. Surfaces of revolution about the Q-axis, for J G {0.0,0.1,0.2,0.3}. The 
radius for given Q is given by the square- root of the cubic —V(Q). For given 
iT", the motion takes place on the intersection of the corresponding surface with a 
plane of constant X. 

Figure 7. Horizontal projection of the solution for an integration of 1000 seconds. 
Top: solution of the 'exact' equations. Bottom: solution of the 'approximate' 
equations. 

Figure 8. Vertical amplitude of the solution for the first modulation cycle (first 
167 seconds). Top: solution of the 'exact' equations. Bottom: solution of the 
'approximate' equations. 

Figure 9. Top panel: Envelope amplitude of the approximate solution. S = 
a/|o| 2 + \b\ 2 (solid line) and C = |c| (dashed line). Bottom panel: Square of the 
eccentricity (solid line) and angular velocity Q (scaled by 50) of the instantaneous 
ellipse (dashed line). 

Figure 10. Azimuth angle (in degrees) for the 'exact' solution solid line) and 
the 'approximate' solution (8, dashed line). The difference 6 — d is plotted as a 
dotted line. The azimuth resulting from integration of ( [4.1 8| ) (not plotted) is 
indistinguishable from the values 9 of the approximate solution. 
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Figure 1. Schematic diagram of the elastic pendulum, or swinging spring. Cartesian 
coordinates centered at the position of equilibrium are used. 
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Figure 2. C is the plane of critical points, A = = B. The vertical axis is R = 
\/\A\ 2 + \B\ 2 . The vertical plane contains heteroclinic semi-ellipses passing from c 

tO — Cq. 
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Figure 3. Phase portrait in the Z-plane for J = 0. The motion is confined within 
the circle \Z\ 2 = |7V. The segment of the imaginary axis within this circle is the 
homoclinic orbit. 
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Figure 4. Plots of Q versus Q for J = and J = 0.25, for a range of values 
E e {-0.0635, -0.0529, -0.0423, -0.0317, -0.0212, -0.0106, 0}. 
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TRICORN SURFACE 




Figure 5. Tricorn surface, upon which motion takes place when H = 0. The 
coordinates are 3 , Q, Q. The motion takes place on the intersections of this surface 
with a plane of constant J (such planes are indicated by the stripes). This surface 
has three singular points. The homoclinic point is marked H.P. 



D. D. Holm & P. Lynch 



Precession of the Swinging Spring 



24 



THREE-WAVE SURFACE, J=0.0 THREE-WAVE SURFACE, J=0.1 




THREE-WAVE SURFACE, J=0.2 THREE-WAVE SURFACE, J=0.3 




Figure 6. Surfaces of revolution about the Q-axis, for J e {0.0,0.1,0.2,0.3}. The 
radius for given Q is given by the square- root of the cubic — V(Q). For given J , the 
motion takes place on the intersection of the corresponding surface with a plane of 
constant X. 




Figure 7. Horizontal projection of the solution for an integration of 1000 seconds. 
Top: solution of the 'exact' equations. Bottom: solution of the 'approximate' 
equations. 
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Figure 8. Vertical amplitude of the solution for the first modulation cycle (first 167 
seconds). Top: solution of the 'exact' equations. Bottom: solution of the 
'approximate' equations. 
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Figure 9. Top panel: Envelope amplitude of the approximate solution. 
S = a/H 2 + \b\ 2 (solid line) and C — \c\ (dashed line). Bottom panel: Square of the 
eccentricity (solid line) and angular velocity Q (scaled by 50) of the instantaneous 
ellipse (dashed line). 
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Figure 10. Azimuth angle (in degrees) for the 'exact' solution solid line) and the 
'approximate' solution (9, dashed line). The difference 9 — ■d is plotted as a dotted 
line. The azimuth resulting from integration of ( |4.18| ) (not plotted) is 
indistinguishable from the values 9 of the approximate solution. 



